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ABSTRACT 


In tills report, we estimate by means of an optimum averaging technique, the 
orbital tracking errors for polar orbiting NASA applications satellites. Accu- 
rate long range estimations of orbital tracking errors are important to users 
in the field as an apriorl condition for orbit updating. The approach we use is 
also called the method of precise conversion of elements, and utilizes results 
of a two week epliemeris produced by a double precision Cowell numerical in- 
tegration calculation, starting with a nominal set of initial orbital elements. 

The gravitational field representation in the numerical integration is complete 
through both order and degree twenty. , Included in this model are the effects 
of solar radiation pressure, lunisolar perturbations, and atmospheric drag. 
Periodic values of inertial Cartesian coordinates of tire satellite state are taken 
from the first week of the epliemeris and used as observational data in an 
orbit determination calculation. Tills data is then used to find the boundary 
conditions or epoch vectors, through a least square processing, of die equations 
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of motion representing analytic solutions >i two separate differential equations 
of un orbit; specifically, those of VintI and Brouwer-Lyddane. These latter con- 
tain only a gravitational field representation through the fourth and fifth zonal 
harmonics respectively, and are well adapted to applications satellite orbits, 
namely low eccentricity, low semi-major axis, and with little or no air drag. 

Both the Vinti and Brouwer-Lyddane methods are then used separately to pre- 
dict an orbit over a period of eleven to thirteen days beyond epoch, and are com- 
pared with the Cowell solution during this time by calculating the in-track errors, 
using Cowell as a standard. The Vinti and Brouwer-Lyddane methods give from 
eight to thirteen kilometers in-track error respectively for the first week of 
free propagation, and approximately twenty-two to twenty-seven kilometers for 
the in-track component after thirteen days. 
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LONG RANGE ORBITAL ERROR 
ESTIMATION FOR APPLICATIONS SATELLITES 

I. INTRODUCTION 

One of the more important tasks associated with the applications satellite sup- 
port systems is pre-flight orbital error estimation. Accurate long term orbit 
prediction is useful as an apriori condition at the field sites for updating the 
orbit. Crucial to the results of such an analysis are not only the particulars 
of the estimation algorithm, but also the physical model used to describe the 
problem under investigation. IIow meaningful these results are will depend to 
a great extent upon how fundamental or comprehensive the system is, as well 
as the accuracy of the estimated initial conditions. A practical test is to make 
a direct comparison with real data from missions of a similar nature, being 
careful to consider the essential differences between the systems. 

In this report, we shall employ a method of optimum orbital averaging or 
precise conversion of orbital elements to study the long range accuracy poten- 
tial of polar orbiting applications satellites such as Landsat and TIROS-N. 

The essence of this approach involves the determination of the boundary con- 
ditions of one set of differential equations of motion, in tiffs case, those 
describing die orbital motion of an artificial earth satellite, by adjusting the 
initial conditions or constants of integration in a least square sense, with the 


use of 'data' generated by another set of differential equations of motion. In 
theory, so that the model might describe a reasonable approximation to the 
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expected behavior of a system sueli as TIROS-N, one needs at the very least, to 
include a data sample commensurate with real cases. For example, one obser- 
vational pair of radio direction cosine data at a specific time Is insufficient to 
determine an orbit, since an indefinite number of ellipses corresponding to this 
piece of information can pass through this point. On the other hand, one Car- 
tesian coordinate and velocity set will completely determine the spacecraft 
trajectory. As a result, it is necessary to process several distinct direction 
cosine pairs in order to fix the orbital elements of the satellite. However, the 
corresponding arc of observations will virtually never exhibit the quality and 
degree of uniformity possible by Interrogating a numerical integration program. 
Likewise, the kinds of uncertainties such as noise, bias, and forces of nature, 
that are inherent in real observational data, cannot be fully and precisely repro- 
duced In a simulation. For the ease in which one wishes to transfer orbital ele- 
ments from a theory that was fitted to real data to any other, the method em- 
ployed in this article then becomes a powerful method of calculation. Uniformly 
spaced, high density data, of the type provided by complete sets of Cartesian 
coordinates would offer one of the best possibilities for exchange of information 
from one set of boundary conditions to those of another orbit determination 
algorithm. 
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In tills paper, we fit analytic calculations to a uniform span of data, spaced 
at hourly intervals, covering periods of ono week and ten days respectively. 
The analytic methods used here are those of Brouwer with modifications for 
small eccentricity and inclination by Lyddane (1), and that of Vinti (2), 

The fitted analytic methods are then used to calculate the orbit for a period of 
two weeks from its epoch, and a comparison with the numerical Integration 
is made by calculating (lie tracking error along tin direction oi motion up to 
thirteen days, 
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II. STATEMENT OF TIIE PROBLEM 


The Satellite Aided Search & Rescue Demonstration will be flown on a TIROS-N 
spacecraft. Di this demonstration, field sites with only minicomputer capabil- 
ity may, under certain adverse contingencies, be required to predict orbits 
for relatively long periods of time. This in turn will directly affect the 
"distress beacon" position location accuracy. In the following calculation we 
consider a nominal TIROS-N orbit. The same general conclusions however, 
apply to other applications spacecraft such as NIMBUS, GEOS, Landsat, etc. 
These nominal Initial conditions for the TIROS-N satellite are as follows: 
a = semi-major axis = 720088G. 36 meters 
e = eccentricity «= zero 
i = inclination = 98. 70 degrees 
M = mean anomaly = zero degrees 
co = argument of perigee - zero degrees 

= longitude of right ascending node = zero degrees 
The epoch time is taken to be January 1, 1977. at zero hours, zero minutes, and 
zero seconds, and the trajectory time span covers twenty-one days. The effec- 
tive cross -sectional diameter and the mass of the spacecraft is taken to be 
3. 489 meters and 711. 687 kilograms respectively. The reflectivity constant 
is 1. 2, while die area for radiation pressure is 9. 560 square meters, and the 
solar pressure constant 4. 50 x 10 “ 3 (kilograms/(second) 2 )/fcilometer. In- 
cluded in the force model of the numerical integration method is the effect of 
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atmospheric drag. The corresponding atmospheric density model is that of 
IIarri8-l 5 riester, 1034, with a profile range between one hundred and one 
thousand kilometers. Lunisolar perturbations are also modelled. The gnivi- 
tational field coefficients are those of the Goddard Earth Modol-I (GEM -I), (3). 

The integration phase was performed using a twelfth order Co we 11- Adams system 
of fixed stepsize, 24 seconds hr duration, since it was anticipated that during 
the compare interval with the analytic theories, the Integration method might 
degrade appreciably as opposed to effects of the gravitational field variations. 

The base inertial system employed In all calculations is that of the mean 
equator and equinox of 1950. 

The method of compare used is described as follows: R N , V N , and L N - 

R n x V N , are the position, velocity and total orbital angular momentum 
vectors of the spacecraft as functions of time determined by the numerical 
integration, R A is the corresponding satellite position vector at the same 
time, and determined by the analytic theory, while the total position error of the 
analytic theory at the given time is the magnitude of die vector difference, 
dR - R N - R a . Then with respect to the numerically integrated orbit plane, 
the instantaneous components of this total position error, along the radial 
direction f\, perpendicular to the orbit plane and along the direction of 

A 

the velocity vector l! , are given by, 

h = fi • AR = Ry_ . AR 
I r n t 


G 


■ ai^= Jx ■ 

I htl 

£ *4 • A R = V N . & R 
respectively. ^ 


The convergence criterion for the orbit differential correction least squares 
solution process is as follows: For a set of m observation residuals, the 
current weighted root mean square value of the observation residual column 

|h |1 £ y [J 

vector at the i iteration is . — —11 — . The corresponding pre dicted weight- 

in 

ed root mean square is H Ay, - A, A *h., II , w here a i s the coefficient 

ni 1 

matrix of the equations of condition, and A x j + 1 is the (1 + 1st) correction 
to the state vector. From this, we have that when, 


Ay, II II Ay, - A, Ax h| il 
in m 

II Ay, II 
m 


<e, 


original rage is, 
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where e 3 UT 4 , the differential correction has converged. 


For the first orbit fitting calculation, both analytic methods use only a gravita- 
tional field representation derived from GEM-1, with Brouwer-Lyddane retain- 
ing those coefficients through the fifth zonal harmonic, and Vinti .hrough the 
fourth zonal harmonic. Jn the second case, both analytic methods employ the 
Goddard Earth Model-5 (GEM-5), (4), representation during the correction 
process, through the same zonal harmonics specified in the first calculation. 
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in. RESULTS AND CONCLUSIONS 

In Tables 1 and 2 wo list tho results of post differential correction comparisons 
of in-track calculations of twenty-four hour Intervals, for both the Brouwer- Lyddano 
and Vlnti analytic orbit calculation methods. Table 1 represents the ease in 
which the analytic theories utilize tho GEM-1 gravitational field representation 
and Tab to 2 gives results Involving the use of the GEM-5 field in tho analytic 
theories. In both tables, time is measured in days, and the in-track error is 
measured in kilometers. From Table 1, the comparison over the period of eleven 
days shows an in-track erro? growth of approximately 2.3 kilometers per daj. 
After one week, this error is approximately 13 kilometers. With use of the 
GEM-5 field the results are somewhat improved (Table 2). The in-track error 
growth for Brouwor-Lyddanc Is approximately 2.1 kilometers per day and for 
Vinti, 1.7 kilometers per day. After one week, Brouwer-Lyddam shows an 11 
kilometer error, and Vinti. 3.5 kilometers. After thirteen days, the in-track 
errors are approximately 27 kilometers and 21.5 kilometers for Brouwer- 
Lyddane and Vinti respectively. Both analytic methods appear to give consistent 


results. 
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Table 1 

Estimated In Track Error For TIROS-N Using 
The GEM-1 Field 


■ 

Time 

(Days) 

Brouwer- Lyddano 
In-Track Error (kilometers) 

Vintl In-Track ' 

Error (kilometers) 

1 

2. 3806 

2. 8083 

2 

3. 74 G1 

3.4503 

3 

5. 4258 

5.2077 

4 

7.5902 

7. 7914 

5 

9.7219 

9.9300 

0 

11.3430 

11. 0218 

7 

13.4270 

12.8735 

8 

10. 4205 

10.4257 

9 

19. 0oj7 

20. 0090 

10 

22.3330 

22.3113 

11 

24. 7967 

- 

24.0057 
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Table 2 

Estimated In Track Error For TEROS-N 
Using The GEM-5 Field 
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In addition, since few noneonservative forces are expected to greatly perturb the 
orbit, ua in-track error growth of approximately 1.7 to 2.3 kilometers per day 
for an orbit of the TIROS-N class, would appear to be a reasonable a priori esti- 
mate. Furthermore, since this growth appears almost linear over eleven to 
thirteen days, one might expect a rapid convergence during subsequent orbit 
improvement calculations. 

As a final consideration in orbit error estimation, let us consider two important 
error sources. If we compare the in-track behavior of a TIROS-N type satellite, 
with no nonconservative perturbations present, using a Brouwer orbit calculation 
with and without the second-order short period terms In the semi-major axis, 
we will get an in-track error growth of approximately 320 meters per day. On 
the other hand, if we assume an uncertainty of one part in 10 6 for u, the product 
of the planetary mass and the gravitational constant, we obtain a value for this 
error growth of slightly more than 300 meters per day. Lyddane and Cohen (5), 
have shown that failure to initially adjust the semi major axis in the second order 
short period terms may immediately show as a serious secular like in-track 
error manifest through the mean anomaly for which 320 meters per day would 
represent a lower bound. This, together with a built in 'boundary uncertainty' 
for the parameter a, might seem to account for the secular in-track error growth 
pattern of Tables I and II. 
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